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The problem of spin-charge separation is analyzed numerically in the metallic phase of the one-band 
Hubbard model in one dimension by studying the behavior of the single-particle Green's function 
and of the spin and charge susceptibilities. We first analyze the Quantum-Monte Carlo data for 
the imaginary-time Green's function within the Maximum Entropy method in order to obtain the 
spectral function at real frequencies. For some values of the momentum sufficiently away from the 
Fermi surface two separate peaks are found, which can be identified as charge and spin excitations. 
In order to improve our accuracy and to be able to extend our study to a larger portion of the 
Brillouin zone, we also fit our data with the imaginary-time Green's function obtained from the 
Luttinger-model solution with two different velocities as fitting parameters. The excitation energies 
associated with these velocities turn out to agree, in a broad range of momenta, with the ones 
calculated from the charge and spin susceptibilities. This allows us to identify these single-particle 
excitations as due to a separation of spin and charge. Remarkably, the range of momenta where 
spin-charge separation is seen extends well beyond the region of linear dispersion about the Fermi 
surface. We finally discuss a possible extension of our method to detect spin-charge separation 
numerically in two dimensions. 

PACS numbers : 71.10.Pm, 71.15.-m, 71.10.Fd, 71.10.Hf 



I. INTRODUCTION 

One-dimensional (ID) interacting fermion systems 
show a number of anomalous properties which cannot 
be understood in the framework of the Fermi-liquid the- 
ory of normal metals. In particular, their momentum 
distribution and density of states are in sharp contrast 
with Fermi-liquid theory for energies and momenta close 
to the Fermi surface. In general, ID systems can be de- 
scribed by an effective low-energy theory based on the 
exactly solvable Luttingcr model with suitably renormal- 
ized parameters and are thus referred to as Luttinger liq- 
uids (LL)JjT3 One of the most striking features of the 
Luttinger model is the complete separation of spin and 
charge degrees of freedom which manifests itself in the 
splitting of the single-particle spectral function in two 
peaks corresponding to sjpttuand charge excitations pro- 
pagating independently .eTou Another important charac- 
teristic of the Luttinger model is the presence of power- 
law behavior with interaction-dependent exponents for 
various correlation functions. 

Beside its application to 1-D systems, LL theory has re- 
ceived particular attention in the past years in the frame- 
work of the theory of high-T c superconductors. The nor- 
mal phase of the high-T c Cu02 planes shows in fact a 
number of anomalous properties which can be possibly 
understood, if one assumes that the CjuCta planes are in 
a kind of two-dimensional LL stateJ3~t!l In particular, 
it has been suggested that spin-charge separation could 
be present also in the Cu02 planes and that it plays an 
essential role in thej-way particles are allowed to tunnel 
between the planes.t2l 



Numerical methods have been proven to be crucial for 
the theoretical understanding of models describing the 
Cu02 planes, since electron correlation is rather strong in 
these systems and perturbative methods are necessarily 
limited. Spin-charge separation is predicted exactly for 
the Luttinger model : an ideal exactly-solvable model. 
It is thus important to test numerically to what extent 
spin-charge separation can occur in a one-dimensional 
physical model. Moreover, in order to prove the theo- 
ries mentioned above, it would be important to check 
whether some two-dimensional models exist, which dis- 
play spin-charge separation. Recently, there have been 
several attempts to detect spin-charge separation in one- 
but also in two-dimensional models. In the U — oo 1- 
D Hubbard modelliij spin-charge separation occurs in a 
natural way at all energies (and not only at low energies 
like expected in a iL) due to the exact factorization p£ 
the wave function.li-3 In a numerical work Jagla et al.t3 
have observed the propagation in real time of a single- 
electron wave packet created at a time t = in a 1-D 
Hubbard model. This wave packet splits up in two exci- 
tations propagating with different velocities that can be 
associated with charge and spin. In a work by W. O. 
Puttika and collaboratorsE^I the possibility of spin-charge 
separation in the 2D t — J model has been signaled by 
the presence of two distinct characteristic wave vectors 
for the spin and charge degrees of freedom. Exact di- 
agonalization of the 1-D t — J modelO have evidenced 
the presence of two peaks in the single-particle spectral 
function whose positions scale with t and J, respectively, 
and have thus been identified with charge and spin ex- 
citations. In another studyE3 two peaks have been de- 
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tected in the single-particle spectral function of a ID 
t — J model with corresponding peaks in the charge and 
spin susceptibilities. These two peaks can be seen, how- 
ever, only for the momenta which are immediately next 
to the Fermi momentum and thus they cannot be asso- 
ciated with a dispepsive spinon and holon band. Finally, 
Kim and coworkers^ have detected two dispersive bands 
for k < kp in a 1-D t — J model close to half filling. 
This is not surprising, since spin-charge separation is in 
fact quite natural to expect when one hole is added be- 
yond half filling. The added hole decomposes indeed in 
a spinless hole and a spin-misalignment which propagate 
with different velocities .E2l Nevertheless, the interesting 
result of that work is that the photocmission spectrum 
of SrCu02, also showing two dispersive bands, is remark- 
ably well reproduced by the numerical results. 

In this work, we present a systematic Quantum- 
Monte-Carlo study of spin-charge separation away from 
half filling, where Luttinger-liquid theory is expected to 
hold, in the whole Brillouin zone (BZ). The nontrivial 
prediction of Luttinger-liquid theory is, in fact, that spin- 
charge separation occurs in the metallic phase, where the 
band dispersion is linear. Spin-charge separation at half 
filling, as studied in the model of Ref. |l9|, is, in our opin- 
ion, of a different nature, since in the insulating phase the 
holon dispersion is quadratic instead of linear. Of course, 
in the case of Ref. USA it was necessary to remain in the 
insulating phase, since the physical system in study was 
half-filled. 

For some values of the momentum k we are able to see 
two peaks in the single-particle spectral function which 
correspond to the spin and charge excitations. However, 
due to the limited resolution of the Maximum Entropy 
method , it is pot possible to resolve the two peaks in 
most of the BZ.E3 For this reason, in the rest of the BZ we 
work with the imaginary-time Green's function Q(k,r) 
which is obtained directly from Quantum-Monte-Carlo 
data without the need of analytic continuation. This has 
the advantage that one does not need to introduce a fur- 
ther source of error produced by the analytic continua- 
tion to real frequencies. Specifically, we perform a nonlin- 
ear x 2 fit of Q(k, t) bv using the solution of the Luttinger 
model (k,rp with two velocities Ui, m, and a 



In Sec. Ill, we discuss and show the results of our fit of 



normalization constant c as fitting parameters. Ell Our fit 
yields a finite value of the difference V2 — v-y larger than 
the statistical error in a large portion of the Brillouin 
zone. Moreover, the fitted values of the corresponding 
excitation energies v\{k — kp) and V2(k — kp) coincide, 
within the statistical error, with the spin and charge ex- 
citations, respectively, calculated independently via the 
associated susceptibilities. It is remarkable that this be- 
havior extends well beyond the region of linear dispersion 
around kp where Luttinger liquid behavior is expected. 

Our paper is organized as follows. In Sec. [ill we intro- 
duce the model, and we show the results of the Quantum- 
Monte-Carlo simulation and analytic continuation to real 
frequencies by means of the Maximum Entropy method. 



the imaginary-time Green's function with the result from 
the Luttinger model. Finally we draw our conclusions in 
Sec. ~ 
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II. QUANTUM-MONTE-CARLO SIMULATION 

We consider the lD-Hubbard model with periodic 
boundary conditions described by the following Hamil- 
tonian: 
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(1) 



where c^ a are annihilation (creation) operators for an 

electron at site i with spin er and m a = c\ a c ia . The 
energy scale t of the model will be set to unity in the rest 
of the paper. 

The simulations were carried out with-. Jthe grand- 
canonical Quantum-Monte-Carlo methodE3'E3 on a 64- 
site lattice with inverse temperature jr^f = P — 20, 
Hubbard repulsion U = 4 and an electron density of 
(n) « 0.75. The simulations yield the one- and two- 
particle Green's functions at discrete imaginary times 
t with < r < (3. We used a discretization of 
the imaginary-time axis At = 0.0625 The spectra 
(one-particle photoemission spectrum, charge- and spin- 
susceptibilities) were then obtained by analytically con- 
tinuing the imaginary-time results to real frequencies by 
means of the Maximum-Entropy methodoCJ'EI. 




FIG. 1. Density plot of the charge susceptibility Xp(9> w ) 
as obtained by the analytic continuation of the Quan- 
tum-Monte-Carlo charge-charge correlation function with the 
Maximum-Entropy method. The grayscale corresponds to the 
value of Xp(<l> w ) (darker regions correspond to larger values of 
X P (q,u>)) and the dots with errorbars show the peak position 
with their uncertainty. The linear fit (straight line) for small 
q yields the charge velocity v p as indicated in the upper left 
corner. 

Figures |l| and || show a density plot of the charge- 
and spin susceptibilities Xp an d Xai respectively. The 
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FIG. 2. Density plot of the spin susceptibility x<t{Q> u ) with 
the same conventions as Fig. |l| The linear fit (straight line) 
yields the spin velocity v a . 



grayscale gives a measure for the value of Xp/o-fe^) as 
a function of momentum transfer q and excitation en- 
ergy uj. The dispersion relation for spin- and charge- 
excitations is defined by the maxima of Xp/a which are 
indicated by dots with errorbars in the figure. A linear fit 
of these maxima near q = yields the spin- and charge- 
velocities v„ = 1.170 ± 0.074 and v p = 2.050 ± 0.093 
which agree very weU— (within the statistical error) with 
Bethe-Ansatz resultsc3 for the infinite lattice and zero 
temperature. 

However, it is not sufficient to have two different ve- 
locities (or, equivalently, energy dispersions) for the two- 
particle spin and charge modes in order to conclude that 
the system shows spin-charge separation. In fact, in a 
Fermi liquid there are spin and charge excitations that 
originate froin. collective modes and do not destroy the 
quasiparticleEj. The quasiparticles thus remain well- 
defined and do not split into a charge and a spin excita- 
tion as it occurs in a Luttinger liquid. On the other hand, 
in a Luttinger liquid (or in spin-charge separated system 
in general) a particle injected at a certain point x decays 
into a spinon and a holon propagating with different ve- 
locities. The separation of the two excitations could then 
be detected by means of a "diagnostic operator" measur- 
ing the time dependence of spin and charge at a given 
point y far away from x. In the case of spin-charge sepa- 
ration, this diagnostic operator would then measure two 
different passing times for the charge and spin perturba- 
tions of the injected particle. True spin-charge separa- 
tion in the sense of the Luttinger model should be thus 
identified with different energy dispersions in the spin 
and charge susceptibilities associated with corresponding 
low-lying excitations in the single-particle spectrum .t3-. 

In Fig. H we plot this single-particle spectrum A(k, w)E2l 
in the whole Brillouin zone. Close to the Fermi momen- 
tum the band dispersion is approximately linear, which 
justifies the mapping to the Luttinger model. However, 
the spectrum becomes broader when going away from 
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FIG. 3. Density plot of the single-particle photoemission 
spectrum A{k,u>) with the same conventions as Fig. jjj. It 
is seen that the dispersion around the Fermi energy (dotted 
line) is linear over a broad momentum range thus justify- 
ing our Luttinger-hquid ansatz for the single-particle Green's 
function. 



the Fermi surface. This phenomenon has two reasons: 
First, the resolution of the maximum-entropy method 
gets worse at higher energies, due to the exponential ker- 
nel in the spectral theorem, and second, according to Lut- 
tinger liquid theory, the single peak starts to split into 
two peaks representing the spin- and charge-excitations 
propagating with different velocities. However, for k very 
close to kp these two peaks, which should be separated 
by an energy («2 — vi)(k — kp), are still too close together 
for the Maximum Entropy method to distinguish them. 
On the other hand, at larger values of (k — kp) the ex- 
citation energies are too high and the maximum entropy 
method becomes less reliable as explained above. In both 
these cases the two peaks merge into a single broader 
peak and spin-charge separation is not detectable. There 
are, however, some favorable intermediate fc-points where 
spin-charge separation is directly detectable in the single- 
particle spectral function. In figure f| we show the spec- 
tral function for one of these favorable points. Here k 
is neither too close nor too far from the Fermi surface 
and the maximum-entropy method (without using any 
prior knowledge) yields two well-separated peaks. Their 
positions are consistent with the spin and charge exci- 
tation energies (indicated by two dots with horizontal 
errorbars) calculated independently from the spin- and 
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FIG. 4. Single-particle photoemission spectrum A{k,u>) 
(in arbitrary units) for k — kF = —4.5 j^. The dots 
with horizontal errorbars indicate the position of spin- and 
charge-excitations calculated by u p / a = (k — kF)v p / a with 
Vp/a obtained from Figs. |l| and |j| For this fc-point close to 
the Fermi momentum the Maximum- Entropy method is able 
to resolve two separate peaks in the spectral function which 
can be identified as the spinon and holon excitation, respec- 
tively. 

III. FIT OF THE IMAGINARY-TIME GREEN'S 
FUNCTION 

In order to carry out a systematic study of spin-charge 
separation it is important to detect spin and charge ex- 
citations over the whole BZ, or at least in an extended 
region around the Fermi surface. However, due to the ad- 
ditional rather large error introduced by the Maximum 
Entropy analytic continuation method to the Quantum- 
Monte-Carlo data, this turns out to be very difficult for 
many k points, as we have discussed above. For this rea- 
son, we work directly with the data for the imaginary- 
time Green's function C7(fc,r). In the asymptotic limit 
(t > 1) and close to the Fermi surface (+fcj?) this func- 
tion should approach the Green's function of the Lut- 
tinger model for right-moving fermions, i. e. 

CX*, ^ r)^Jdx e- ik *Gl™] Kp {x , T ), (2) 

with 



y/VlT + lXy/v 2 T + IX{X 2 + V 2 T 2 ) ~ ( K P + 1 / K P " 2 )/ S ' 

where c is a normalization constantil and kp the Fermi 
momentum. Therefore, in order to identify the spin and 
charge excitations directly in the Green's function, we 
carry out a nonlinear \ 2 fit °f our data for Q{k,r) to 
g\ LM ) (k.r). The fit parameters are the two velocities 

v\ and V2 , and the normalization constant c Ej Due to the 
statistical error in the Quantum-Monte-Carlo data, we 



get statistical errors A^i , At>2 , and Ac for the parameters 
Vi, V2, and c, respectively. The splitting of the single- 
particle mode into two excitations is thus detected when 
the difference between the two velocities is larger than the 
statistical error. Furthermore, in order to make sure that 
the two excitations coincide with the spin and the charge 
modes one has to compare v\ and vi with the velocities v p 
and v a calculated independently via the susceptibilities. 
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FIG. 5. Logarithmic plot of the imaginary-time Green's 
function Q(k,r) vs r with k — kF + 7r9/64, as obtained from 
Quantum-Monte-Carlo data. 

However, in order to carry out this fit one should not 
use the data from the whole interval 0.0 < r < (3 for 
the following reasons. First of all the Hubbard model 
and the Luttinger model Green's function should coincide 
only asymptotically. For this reason, we choose r > 1.0. 
Moreover, the interval /3/2 < r < (3 is equivalent to the 
one — (3/2 < t < Oso that we can omit the former. In ad- 
dition, as can be seen in Fig. |5| the log of the imaginary- 
time Green's is quite sharply defined up to r ~ 5.0. For 
t > 5.0, large (relative) errors start to develop due to the 
small value of the Green's function in these points. For 
this reason, we choose to carry out our fit only for the 
data in the interval 1.0 < r < 5.0 in order to select the 
less "noisy" data. In order to check that our results do 
not depend on this choice, we also carry out the fit for 
the data in the interval 1.0 < r < (3/2 = 10.0 (Fig. @). 
This turns out to be quite similar to the first one (Fig. 
^) except for larger statistical errors. In Fig. || we show 
the result of our fit with the T = Green's function [ (|^) 
with (||)] for several values of q = (k — kp)- The ver- 
tical black lines show the value of the spinon and holon 
excitation energies £i{q) = vi(q) q and £2(9) = V2(q) <?, 
respectively, obtaiaed from the fit with the single-particle 
Green's function.ca As one can see, we obtain a clear sep- 
aration of the two modes for almost all the q points. In 
addition, the velocities are slightly g-dependent as ex- 
pected from a curved band. To check that these modes 
correspond to spin and charge degrees of freedom, we plot 
in the same figure the dispersions calculated from the 
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peaks of the corresponding susceptibilities. The width of 
the gray regions indicate the peak positions within their 
uncertainty. As one can see, the dispersions obtained in 
the two ways coincide within the statistical error. We 
find it remarkable that, even at fc-points far from fcj?, 
the fit with the Luttinger-liquid Green's function is in 
agreement with the two-particle response, although the 
dispersion is no longer linear. It thus seems that spin and 
charge separation survives even at higher energies. 




FIG. 6. Spin and charge dispersions ei and £2 vs 
q — k — Hf (errorbars without and with central dot, re- 
spectively) as obtained from the x 2 fit of the Quan- 
tum-Monte-Carlo data for the imaginary-time Green's func- 
tion with the Luttinger liquid Green's function [ with ([])]. 
The Luttinger liquid Green's function is taken at zero tem- 
perature and with correlation exponent K p — 1. The fit is 
carried out for the data in the time interval 1.0 < r < 5.0. 
For comparison, we also show the dispersions obtained from 
the peak positions (with corresponding uncertainty) of the 
spin (dark gray) and of the charge (light gray) dispersions 
(Cfr. Figs. |and|). 

In Fig. ^ we used the simplest form of (||), namely, 
the one with K p = 1.0. This is not, in principle, the 
correct value of the correlation exponent K p when U =/= 0. 
Actually, one could use K p as a further parameter to fit 
the data or, altewiatively, use the result from the Bcthe- 
Ansatz solutionEl It turns out, however, that an attempt 
to fit K p yields an error of the order of 0.5, which means 
that Kp cannot be determined by our fit. It also turns 
out that the result of the fit does not depend crucially 
on the value of K p we are using. Indeed, as one can see 
from Fig. ||, where we show the results of the fit obtained 
with the Bcthc-Ansatz value K p — 0.7 the results of the 
fit do not differ appreciably from the ones in Fig. [|. For 
this reason, the non-interacting value K p — 1.0 can be 
safely used. This is important, because in this way it is 
possible to test the occurrence of spin-charge separation 
even without knowing whether the system has anomalous 
scaling {K p 7^ 1) or not. This could be useful in cases 
where the anomalous exponent may be not known a priori 
and may be difficult to evaluate. In this case, one can 
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FIG. 7. Same as Fig. q except that the fit is carried out 
for the data in the time interval 1.0 < r < 10.0. 
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FIG. 8. Same as Fig. ^ except that the fit is carried out 
with the Luttinger liquid Green's function [ (|^) with (^)] with 
correlation exponent K p — 0.7. 



assume a form of the fitting function without spectral 
anomaly (i. e., K p = 1), but simply with a branch cut 
due to spin-charge separation. This could be useful, for 
example, to test spin-charge separation in 2-D as we shall 
discuss below. 

Since the Quantum-Monte-Carlo simulations are car- 
ried out at finite temperature (3 = l/(fcsT) = 20, we 
also perform our fit with the Luttinger model Green's 
function at the same temperature Q^fti k r ' @ = ^0), 
which has the same form as (||) with (j|) replaced with 



(4) 
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with 



UtK (x,T;f3) 



(5) 
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cos(f) 



K+l/K 



-i arg(tanh x+i tan t) 



In Fig. ^, we show the fit performed with the more 
complicated finite-temperature (/? = 20) Green's func- 
tion (||). As one can see, the results are not appreciably 
different from the ones of Fig. ^, which means that the 
temperature of our simulation is low enough and we can 
safely fit our results with the T = Green's function. 
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FIG. 9. Same as Fig. |^ except that the fit is carried out 
with the Luttinger liquid Green's function [ (12) with (El)] with 



finite temperature (/3 = 20) and correlation exponent K p 



1. 



Finally, to check that the two different velocities v\ and 
i>2 obtained are not an artifact of our fit, we carry out a fit 
of the non-interacting Green's function (k,r; /3) with 
the function (||) assuming artificially the same statistical 
errors as the ones obtained in the Quantum-Monte-Carlo 
simulation. As one can see in Fig. [h], in this case the 
spin and charge velocities obtained are equal within the 
statistical error, as it should be. 
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FIG. 10. Fit of the U — imaginary-time Green's function 
with Eq. (^). The straight line shows co — vpq where vp is 
the Fermi velocity of the U — system. 



Another motivation of this work was to test a "diagnos- 
tic operator" that can be applied to detect numerically 
the occurrence of spin and charge separation in a many- 
body system from Quantum-Monte-Carlo data. If one 
uses exact diagonalization, where the spectral function 
of the system can be evaluated directly, it is of course 
not necessary to fit the imaginary-time Green's function. 
However, we believe that the systems that can be stud- 
ied by exact diagonalization (10-16 sites) are too small 
to allow for a systematic study of spin-charge separa- 
tion (except for a very high value of the momentum, like 
in Ref. |l5). The Fourier transform (in momentum and 
imaginary-frequency space) of the spin-charge separated 
Green's function Eq. (||) with K p — 1 reads 



1 



y/ito- e\{k)y/iuj - e 2 (k) 



(6) 



where £i(fc) = t?i(fc — kp) and £2(fc) = v 2 {k — kp) rep- 
resent the spin and charge excitations (measured from 
the chemical potential) in which the single-particle exci- 
tation is split. The same form of G^^k s i w ) could be 
expected to hold asymptotically, i.e. for small frequen- 
cies and close to the Fermi surface, in higher dimensions. 
Close to the Fermi surface, one will have a direction- 
dependent dispersion £j(fe) = (k — kp) ■ v(kp) where 
v(kp) is the Fermi velocity of the point at the Fermi sur- 
face kp closest to k. Spin-charge separation would be 
signaled by two different, direction-dependent £i(fc) and 
£2(fe) for a given k. 



IV. CONCLUSIONS 

To summarize, we carried out a test of spin-charge sep- 
aration in the 1-D Hubbard model at finite doping. It 
is in general difficult to resolve the peaks corresponding 
to the spin and charge excitations in the single-particle 
spectral function due to the loss of accuracy which oc- 
curs when analytically continuing the imaginary-time 
Quantum-Monte-Carlo results to real frequencies. For 
some values of the momentum close to the Fermi surface, 
however, we were able to resolve two peaks whose ener- 
gies correspond to the peaks at the same q — k — kp in 
the spin and charge susceptibilities, respectively. 

By fitting the Quantum-Monte-Carlo data for the 
imaginary-time Green's function with the exact solution 
from the Luttinger model with the spin and charge veloci- 
ties as fitting parameters, we have been able to resolve the 
two excitations over the whole Brillouin zone. The two 
excitation energies found in the fit agree, within statis- 
tical error, with the spin and charge excitations, respec- 
tively, identified with the peaks of the spin and charge 
susceptibilities. Remarkably, this occurs also away from 
the region where the band dispersion is linear. We also 
suggested a possible extension of this "diagnostic opera- 
tor" to test a possible occurrence of spin-charge separa- 
tion in two dimensions. 
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